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Abstract. We reconsider parity violation experiments in atomic hydrogen and 
deuterium in the light of existing tests of the Electroweak interactions, and assess 
whether new experiments, using improved experimental techniques, could make useful 
contributions to testing the Standard Model (SM). We find that, if parity experiments 
in hydrogen can be done, they remain highly desirable because there is negligible 
atomic-physics uncertainty and low energy tests of weak neutral current interactions 
are needed to probe for new physics beyond the SM. Of particular interest would be 
a measurement of the nuclear spin independent coupling Cid for the deuteron at a 
combined error (theory -|- experiment) of 0.3%. This would provide a factor of three 
improvement to the precision on sin^ Ow at very low momentum transfer provided 
by heavy atom Atomic Parity Violation (APV) experiments. Also, experiments in 
H and D could provide precise measurements of three other electron-nucleon, weak- 
neutral-current coupling constants: Cip, C2p, and C2D, which have not been accurately 
determined to date. Analysis of a generic APV experiment in deuterium indicates that 
a 0.3% measurement of Cid requires development of a slow (77K) metastable beam of 
« 5 X 10^^D(2S)s~"'^ per hyperfine component. The advent of UV radiation from free 
electron laser (EEL) technology could allow production of such a beam. 
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1. Introduction 

In this report, we answer two questions concerning measurement of atomic parity 
violation (APV) in hydrogen and deuterium. First, given the current level of tests 
of the Standard Model (SM), how well would parity violation need to be measured in 
hydrogen and deuterium to make a useful contribution to this body of data? Second, 
what experimental parameters are needed to achieve the required precision? 

As to the first question, there exists a long list of impressive experimental tests of 
the SM [I], many involving large collaborations in high energy physics. However, as 
discussed in section [21 future progress may, in part, require additional high-precision, 
low-energy measurement s|2j, and in this regime, experiments in hydrogen have the 
advantage that the atomic structure is known to high precision and would not limit 
the theoretical interpretation. 

Another advantage of parity experiments in hydrogen is that all of the electron- 
nucleon, weak-neutral-current couplings can be determined by performing measurements 
at several level crossings in both H and D. This point is discussed in Section [3] where 
we present updated SM predictions for the parity violation mixings at the various 
crossings. A summary of the results of former PNC measurements in hydrogen is 
presented in Section HI and in Sections [5], and [6] we consider our question concerning 
the experimental requirements by evaluating a generic APV experiment in deuterium. 
Finally, we summarize our conclusions in section [71 

2. Theoretical Motivation 

Testing the SM of the electroweak interaction remains an important task of experimental 
physics. This effort has reached the level of precision measurements. For example, the 
weak mixing angle sin^ 9\y has been determined to 0.06% at the Z boson mass. Most of 
the data confirming the SM have been obtained at high energy near the Z resonance. 
The excellent agreement at the Z-pole was obtained by the weighted average of two 
experiments that disagree with one another by 3.4 a. Moreover, there have been recent 
measurements at lower energy, such as the muon g-2 experiment [3], the NuTeV neutrino 
deep inelastic scattering experiment [IJ , and the SLAC E-158 polarized M0ller Scattering 
experiment [5]. In the muon g-factor experiment there is a 3.4 standard deviation 
difference between the experimental result and the theoretical calculation [6] suggesting 
new physics beyond the Standard Model. This situation has motivated additional low 
energy measurements of the electroweak interactions, such as the Qweak [71 [8] and 
DIS-Parity [9l [10] experiments at the Thomas Jefferson National Accelerator Facility 
(JLab). 

A useful way to assess the low energy tests of the electroweak interaction is to 
consider the predicted running of weak mixing angle sin^ 9w- The SM predicts that this 
parameter evolves as a function of momentum transfer [21 [Til [12]. This evolution can be 
checked by performing precision low energy experiments and comparing the results with 
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precise measurements at high energy. In this respect, tests done near zero momentum 
transfer such as APV or low energy electron scattering experiments are ideal. In the 
MS scheme, Erler and Ramsey-Musolf [21 [12] find for the value of the weak mixing 
angle at zero momentum transfer: 

sin^ OwiO) = 0.23867 ± 0.00016. (1) 

In the same scheme the value at the Z-pole is [T]: 

sin^ 9w{Mz) = 0.23122 ± 0.00015, (2) 

so that the zero momentum transfer value is larger by about 3%. Several low energy 
measurements test this prediction. APV (Cs[T3] and Tl|14t [T^) establishes the running 
of weak mixing angle with a 0.9% measurement of sin^ 9w that is one standard deviation 
below the SM prediction. The E-158 M0ller scattering experiment has an error of 
0.6% which is one sigma above the SM. NuTeV, a deep inelastic neutrino scattering 
experiment, reports a 3cr discrepancy with the SM value. 

The Qweak experiment aims to improve the limits on the semi-leptonic weak neutral 
current couplings by measuring the weak charge of the proton Qwiv)- In lowest order 
the weak charge is given by [16]: 

Qwip) = (l-4sin2^H/). (3) 

Since sin^ ^vk is close to 1/4, Qw{p) is suppressed and particularly sensitive to higher 
order corrections and new physics beyond the SM. 

To be more precise, we can use the SM expressions for the electron-hadron 
processes from Ref. [1], including radiative corrections, to derive the relation between 
the uncertainties in the weak charge Q\y{Z,N) and sin^^^ for an atom with nuclear 
charge Z and neutron number A^. We find: 



A(sin'^VK) / ^^^o 1 n^70^^\ 



-0.078 + 1.0785^ ) (4) 



sin Ow \ Z J Qw 

For H, D, and Cs the factors multiplying AQvv/Qiy are -0.078, 1.00, and 1.45, 
respectively. For example, a measurement of Qw{p) with a precision of 4% determines 
sin^ 9w to 0.3%. For the Boulder Cs experiment, if the experimental error (~ 0.35%) 
[T3] and atomic theory error (~ 0.5%) are combined in quadrature, the total uncertainty 
in Qw{Cs) is ~ 0.6% which corresponds to an uncertainty of ~ 0.9% in sin^ based 
on Eq. [D, and this is the dominant contribution to the current APV limit on sin^ 9w- 
For deuterium, Eq. [H indicates that the fractional uncertainty in Qw{D) is equal to the 
fractional uncertainty in sin^ 9w 

The need for precision measurements of semileptonic weak neutral current couplings 
at low momentum transfer obliges a renewed consideration of measurements in hydrogen 
and deuterium where precision is not limited by calculational uncertainties. In this 
respect, both the heavy atom APV experiments and polarized electron scattering 
experiments are susceptible to hadronic or atomic effects. Qweak aims to measure the 
parity violating elastic electron-proton scattering asymmetry at ^ 0.03 GeV^. This 
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process is proportional to the sum of Qwiv) and a proton form factor F^{Q'^,6). The 
form factor contribution must be subtracted from the scattering asymmetry in order to 
extract Qw{v)- In principle, this extrapolation would require a hadronic structure 
calculation, but an extensive series of parity violating electron-proton scattering 
experiments at higher will allow extrapolation to (5^=0 with an uncertainty of about 
2% in the weak charge[2]. The heavy atom APV experiments require sophisticated 
atomic theory calculations in order to extract the weak charge from the experimental 
results. In the case of Cs, the atomic structure uncertainty is currently about 0.5%. To 
go beyond this in heavy atoms probably requires measuring APV for a series of isotopes 
of an atom and forming ratios of the results. These ratios determine Qw and are less 
sensitive to calculational uncertainties. However, this procedure is subject to other 
uncertainties arising from variations in neutron distributions along the chain, which will 
limit the ultimate precision [2J. 

Apart from the potential of a precision measurement of sin^ O-^y and sensitivity 
to new physics, experiments in hydrogen and deuterium provide the opportunity to 
determine all four of the electron-nucleon, weak-neutral-current coupling constants: Cip, 
C2p, Cin, and C2n (13 HE]. Here Cip and Ci„ are the nuclear spin independent couphngs 
for the proton and neutron. (Note that Qwip) = 2Cip [HI [H].) and are 
the nuclear spin dependent couplings. These four parameters are important because 
they are basic properties of the neutron and proton. The constants for the neutron 
can be obtained by comparing results from hydrogen and deuterium experiments. For 
deuterium we use the notation Il9i : 



In terms of the quark coupling constants Ci„, Cid, and C2d (Ref. [T]), the proton 
and deuteron couplings have the values |20l [21], [221 [23] : 



The constants Cip and Cid can be determined from APV experiments with little strong 
interaction uncertainty. By contrast, C2p and C2D involve considerable strong interaction 
uncertainty, as has been discussed by a number of authors[22l [2H [25] . 

In Table [1] we give experimental and theoretical values for these couplings. The 
experimental numbers in column 2 are based on a "model-independent fit" to all neutral 
current parameters which are presented in Table 10.8 of Ref. [T]. The fit allowed for an 
arbitrary electroweak gauge theory. The final column of Table [Tjgives the SM predictions 
based on Eq. [6l Table 10.3 of Ref. [1], and siv? 9w = 0.23122. Theoretical corrections 
from the anapole [26] and strange quark [27] contributions are not included. All of 
the couplings except Cic are suppressed in the Standard Model. For Cip and C2p the 





(5) 




(6) 
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Table 1. Parity violating weak neutral current couplings for H and D. 



Coupling From model-independent fit[lj Standard Model 



Cip 


0.082±0.034 


0.0355 


CiD 


-0.441±0.012 


-0.4587 


C2p 


-0.16 ±0.26 


0.043 


C2D 


-0.32 ±0.4 


0.007 



suppression arises because they are proportional to 1 — 4 sin^ 9w in lowest order. The 
coupling C2D vanishes exactly in lowest order and the small residual value is entirely 
due to higher order corrections, so this parameter is quite sensitive to these corrections. 
Although CiD has been determined (indirectly) to about 3%, the uncertainty in Cip is 
comparable to the size of the SM prediction, and both C2p and C2D are undetermined. 

An attractive goal for a future APV experiment would be a 0.3% measurement 
of the nuclear spin independent coupling for the deuteron Cid- This would determine 
sin^ 9w{0) with a combined error (theoretical ± experimental) of 0.3% which is a factor 
of three improvement on the limit provided by the heavy atom APV experiments. This 
is comparable to the limit on sin^ 9w{0) expected from the Qweak experiment, but 
would provide complementary information on the isospin dependence of the e-nucleon 
electroweak couplings [28] . 

A measurement of the nuclear spin independent coupling in hydrogen at the same 
level of absolute precision would provide a 4% measurement of Cip, which together 
with a measurement of Cid would provide a constraint on the isospin dependence of the 
nuclear-spin-independent weak-neutral-current interaction with one experimental setup. 
For C2p and C2D, measurements with absolute precisions similar to a 0.3% measurement 
of CiD could provide the first determinations of these quantities with an error of 3% 
for C2p and 20% for C2d- Again, the isospin dependence of the nuclear spin dependent 
couplings could be constrained by comparing the two measurements. 

3. Parity mixing in hydrogen and deuterium 

In planning parity experiments in hydrogen, it is natural to take advantage of the 
metastable 2^Si/2 state. Metastable hydrogen atoms live long enough (~ 122 ms) to 
survive the length of a typical laboratory atomic beam apparatus, and because of the 
near degeneracy of the 2^Si/2 and 2^Pi/2 states, level crossings occur at relatively modest 
magnetic fields (~ kG), and this provides a means for separately determining the weak 
interaction coupling constants, since different combinations of constants are resonantly 
enhanced at the various level crossings. In order to determine all four coupling constants 
in Table [H it is necessary to perform experiments at several magnetic field settings in 
both H and D, and combine the results to isolate the separate constants. In this section, 
we review the parity mixing matrix elements at the various crossings in H and D and 
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Figure 1. Zeeman diagram sliowing the 2^Si/2 and 2^Pi/2 levels in hydrogen in 
a static magnetic field|29). The states are labeled using Lamb's nomenclature [30] . 
Subscripts give the values of the good quantum number mp. 



present updated SM predictions for their strengths. 

Fig. [1] shows the relevant energy levels of hydrogen in a static magnetic field of 
400-1400 G. In this region, the Zeeman energy from the coupling to the magnetic field 
is much less than the fine structure splitting so the total electronic angular momentum 
J is still a good quantum number. On the other hand, at fields above ~ 200 G, the 
electron and nuclear spins are decoupled so the natural basis is |JmjImi), where I is the 
nuclear spin. In Fig [1] the 2^Si/2 levels are labeled a (mj=+l/2) and /? (mj=— 1/2), 
while the 2^Pi/2 states are labeled e (mj=-|-l/2) and f (mj=— 1/2) [31]. The subscripts 
are mp = mj + mi. 

Fig. [2] shows blowups of the jSe level crossings in hydrogen and deuterium. The 
selection rule for HpNc is Amp = 0, so only the levels Po and cq in hydrogen are mixed. 
The value of the matrix element is given in Table [21 which indicates that a measurement 
at the Pe crossings in hydrogen is sensitive only to the nuclear spin dependent coupling 
C2p. In deuterium, parity mixing is resonant at two /5e crossings which, as indicated in 
Table [21 are sensitive only to C2d- 

Fig. [3] shows the /3f level crossings in hydrogen and deuterium. Parity mixing is 
resonant at two crossings in hydrogen. Table [2] indicates the matrix elements depend 
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Figure 2. Parity mixings resonant at the /3e crossings in hydrogen (upper part) and 
deuterium (lower part). 



on nearly orthogonal combinations of C\p and C2p- So parity experiments designed to 
use the /3f crossings in hydrogen provide a means to measure both of these constants. 
For the /3f crossings in deuterium, there are three crossings of levels mixed by the weak 
interaction. From Table [2] it is seen that two of these depend on orthogonal combinations 
of C\B and C^d while the crossing of /-1/2 and /3-i/2 is, to good approximation, sensitive 
only to C\D. 



4. Current limits on parity violation in hydrogen 

There have been numerous suggestions for APV experiments in hydrogen. One 
possibility is to observe PNC asymmetries in the single photon decay 2S1/2 IS1/2 
+ 7 of a metastable beam in vacuum [T7]. or in a beam perturbed by applied electric 
and magnetic fields [32l [33] . These ideas have not led to practical experiments in 
hydrogen although they hold promise for experiments in heavy one-electron ions[33]. For 
hydrogen, the most promising approach is to drive optical or microwave transitions 2S1/2 
— *• nSi/2 and look for a pseudoscalar dependence on the transition rate[3 ^ [35 |[36l[37| [38]. 
Microwave (or RF[39j) transitions 2s — 2s' within n=2 are particularly attractive 
because Doppler broadening is negligible, and linewidths can be small. Also, relatively 
high power is readily available and the interaction volume can be large. The possibility 
for narrow linewidths is particularly significant since it allows complete separation of the 
various hyperfine transitions, and, as discussed in Section |5|, a narrow line is required 
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Figure 3. Parity mixings resonant at the j3f crossings in hydrogen (upper part) and 
deuterium (lower part). 



Table 


2. Matrix elements for parity-violating, 


weak-neutral-current couplings for 


hydrog 


en and deuterium at the level crossings. 


In units of iVw w 27r x 0.013 5"^ 


ClD = 


Cip + Cin and C2D = C2p + C2n- 




Matrix Eloiiieiit 


Value 


Standard Model 


Hydrogen 






{eo\V\(3o) 


— 2C2P - 


-0.086 


{h\v%) 


Cip + I.IC2P 


0.083 




Clp — C2p - 


-0.007 


Deuterium 






Mv\Pi/2) 




-0.010 


(e_i/2|^|/3-i/2> 




-0.010 


(/l/2|^|A/2) 


ClD + C2D - 


-0.452 


{f-l/2\V\(3^l,2) 


CiD - 


-0.459 


(/_3/2|l^|/?-3/2) 


ClD — C2D - 


-0.466 
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to optimize the sensitivity to parity violation. 

Starting in the late 1970s, hydrogen parity experiments were performed at the 
University of Michigan [311 EEl iQl iH Hi], The University of Washington [3H1 US], 
and Yale University [35l [39]. The experiments were motivated in part by the near 
degeneracy of the 2^Si/2 and 2^Pi/2 levels in hydrogen which enhances the parity mixing 
and provides an opportunity for detection of PNC effects even though hydrogen does 
not benefit from the enhancement as a function of nuclear charge which motivates 
parity experiments in heavy atoms [46]. 

All of the hydrogen APV experiments were based on fast (~ 500 eV) beams 
of metastable hydrogen atoms and aimed to measure C2p- Typically, the beam was 
prepared in one of the a hyperfine states and a transition a —>■ /3 to one of the (5 
hyperfine states was driven. The transition was observed by selective quenching of the 
resulting /? state atoms and detecting the Lyman-a radiation emitted. 

The transitions were driven by a parity conserving amplitude Apc and a parity 
nonconserving amplitude ApNc- The parity nonconserving amplitude proceeds via an 
intermediate 2Pi/2 state: 

e"^ Vw 

a 2Pi/2 ^ /?, (7) 

where the first step is an El transition driven by microwave electric field e'^, and the 
second step is the 2^Si/2 — 2^Pi/2 coupling by the weak interaction. This is also called 
a "weak induced amplitude". 

Two types of parity conserving amphtudes Apc were used. One of these was an 
Ml amplitude [35]: 

a 13, (8) 

where If is the microwave magnetic field. More commonly a "Stark induced amplitude" 
was used: [361 138] : 

a 22Pi/2 (3, (9) 

here, the first step is a microwave El coupling to an intermediate 2^Pi/2 state driven 
by e"^, and the second step is a 2^Si/2 — 2^Pi/2 coupling by a static electric field E'^ . 

The experiments searched for a parity violating (pseudoscalar) interference term 
between ApNc and Apc in the transition rate, that changed sign as the "handedness" of 
the interaction region was switched. The main differences among the various schemes 
were in the configurations of fields in the interaction regions, and the pseudoscalars of 
interest. Particularly significant is that in the Michigan experiment, Apc and ApNc were 
driven simultaneously in a single microwave cavity, while in the Yale and U. Washington 
experiments, separate oscillating field regions were used for ApNc and Apc. 

These experiments were not able to observe parity violation (for a short review see 
[47]). The best limit was obtained by Fehrenbach [44] who reported: 

C2p = 1.5 ± 1.5 ± 22. (10) 
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Here, the first error is statistical and the second is systematic. From Table [H we see that 
the value expected based on the Standard model is C2p= 0.036. So the statistical error 
is a factor of 40 larger than the expected effect while the systematic error is 600 times 
larger than the expected effect. It is clear that more than an incremental improvement 
is necessary if a useful measurement is to be obtained. Some ideas of how this might be 
accomplished are given in sections [5] and [61 



5. Sensitivity of a generic deuterium APV experiment 

The original hydrogen parity experiments could not have observed PNC. The problems 
were a lack of sufficient sensitivity and large systematic errors. In the next two sections 
we discuss a particular idea for an APV experiment in deuterium aimed at providing a 
0.3% measurement of Cid- The conclusions of this paper do not rest on this particular 
scheme, as there are many other possibilities which would have roughly the same 
sensitivity. Rather, we chose this example after considering some general guidelines for a 
future H/D parity experiment. First, in comparison to the original H APV experiments, 
the signal-to- noise-ratio S/N has to be improved by more than two orders of magnitude. 
Then, to improve handling of systematic effects, it is important to have precise control 
over the alignment of the apparatus and all applied fields. To make this less daunting, 
the apparatus should be simple - a few homogeneous apphed fields is preferred over 
more complicated configurations. Perturbations of the beam due to surfaces (apertures, 
walls, etc.) should be minimized, and all surfaces should be kept as far from the beam 
as possible. The apparatus should be "clean" i.e. surfaces should be carefully prepared 
to minimize stray electric fields, and ultra-high vacuum techniques should be used. The 
beam density should be as low as possible consistent with the needed intensity and 
there should be no charged particles, or field configurations which tend to trap charged 
particles. 

One can increase S/N by developing a more intense metastable beam. In addition, 
there are advantages in both sensitivity, and in reducing systematic errors to using 



a slow beam [381 1^ . In [Appendix B| we discuss the requirements for producing a slow 
(77K) metastable beam with about 5 x 10^^ D(2s)s~^ per hyperfine level. An important 
consideration is that there is a limit on the density of metastable atoms in the beam 



due to quenching by beam-beam collisions. This is discussed in [Appendix A 

A particularly simple configuration for the interaction region would involve driving 
an Ml transition [35] with a plane polarized microwave beam incident at right angles 
to the metastable beam. This would eliminate the need for an additional applied 
electric field to drive Apc. Also, eliminating the microwave cavity (or cavities), has 
some advantages. The standing wave pattern in a cavity leads to spatial variations 
in the fields which can complicate the transition amplitudes [H] and make analyzing 
systematic errors more difficult. In addition, apertures at the entrance or exit to the 
cavity perturb both the atoms and the cavity modes. 

A practical way to realize this is to use a microwave-open-ring- resonator [l9l [50|, [51] 
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X 




Figure 4. Fields in the interaction region. The metastable D beam is traveUing along 
the z axis. The microwave electric e and magnetic b fields lie in the x-z plane. B is a 
static magnetic field. 

in the interaction region. Complete control of the polarization in such a resonator is 
possible using quasi-optical techniques [5^. Such a resonator allows the propagation 
direction of the traveling waves to be switched, providing the ability to reverse the 
microwave propagation vector k. 

As an example, we consider a measurement of Cio at the crossing of the levels 
f-i/2 and (3-1/2, utilizing the configuration of fields shown in Fig. HI The D(2s) beam 
is directed along the z axis, coUinear with a static magnetic field B. The beam of 
microwaves propagates into the page (k = — y) and the frequency is set to drive the 
Ml transition a+1/2 — ^ /5-i/2- The microwave electric polarization vector e lies in the 
x-z plane, and is orthogonal to the microwave magnetic field b. The vector e is a 
combination of fields and [52j . These fields are derived from the same source but 
have a phase difference of 6 and separately controlled amplitudes. The resulting field is 
given by: 

e(t) = e^; cos(z/t) + €2 cos(z/t + 6^). (11) 

Here, we will set 6 = or tt, in which case, the microwave polarization is linear. The 
polarization angle relative to the x axis is determined by the relative strength of 
and Sz- When 6 is changed from to vr, goes to -(f). Other values of 6 give circular 
or elliptical polarization. The amplitudes in the experiment are shown schematically in 
Fig [5l The parity conserving part Apc is an Ml amplitude driven by the component of 
the microwave magnetic field along the x axis (b^;). 

The parity nonconserving amplitude ApNc proceeds via the intermediate level 
f-i/2 and involves the weak interaction HpNc and an ii^l coupling via the microwave 
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Figure 5. Parity conserving Apc and nonconserving ApNc amplitudes. 



electric field component along the x axis (e^;). The analysis of this experiment requires 
consideration of the equations of motion for the three levels a+1/2, P-1/2, and /-1/2 but, 
to sufficient approximation for the present considerations, can be reduced to an effective 
two-level problem by introducing the parameter ja+-^/2 which is the decay rate of the 
level due to the quenching a+1/2 f-1/2 by the microwaves [IT]. The amplitudes 
a for the state and b for the P-1/2 state obey the equations: 

exp[-iK^^/2^_^/^ - u)t]b - lha+,^,a ^^^^ 

Aexp[+i(u;„^^/^^_^/^ - u)t]a 



la 
ib- 



where [531: 



2hA 



-/iobx - iVwCiD(i2esA/3_^/2f_j/2. (13) 

The parameter = 27r(0.013)s^^ gives the magnitude of the weak interaction matrix 
element, and d2=V^eao- The parameter ^i3_i^2f-i/2 defined by: 

1 



A, 



UJ. 



(14) 



UJ 



and 72p = 27r x lO^s ^ is the natural decay rate of the 2p state. 



Here, uj. 

The beam is prepared in the 0+1/2 state so the initial conditions are: 
a(0) = 1 
6(0) = 0. 

For simplicity, we treat the case of exact resonance cUa^^^/j/^-i/- ~ ^ 



(15) 



0, and also choose 



=0, 



/2/-1/2 



the maximum interference between the two terms in Eq. [iniwhich requires 
i.e. that the magnetic field be tuned to the center of the level crossing. If the atom 
spends a time r in the interaction region, the amplitude to be in the state /5_i/2 after 
this time is given by: 

1 -exp(-7„ r/2) 



b{r) 



-iA- 



(16) 



The two terms in A (Eq. [13]) determine the two amplitudes of interest. Apc is 
proportional to the first term while ApNc is proportional to the second. Since the 
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Table 3. Parameters determining the optimized sensitivity (Eq. [22]) 



Symbol 


Meaning 


Assumed for Figs. [6] and [7] 


V 


Total detection efficiency 


0.5 


J 


Metastable flux/hyperfine component 




T 


Counting time 


5 days 


T 


Interaction time 




i% 


H PNC coupling strength (n=2) 


2tt X 0.013 s-^ 


l2p 


2p decay rate 


27r X 10^ s-^ 



electric e and magnetic b fields in the plane wave are orthogonal, = bsin0 and 
ex = e cos 0, so the sign of the interference term changes as — —0. The interference 
term is proportional to the pseudoscalar (e ■ B)(e x k) ■ B, so it also changes sign if the 
microwave propagation vector is reversed, i.e. k — > — k; but it is even under a reversal 
of B. This pseudoscalar can appear in a T-even theory due to the presence of damping 
of the 2p state [371 ill |13] . 

The parameter 7q^^/2 determined in a separate calculation. If the magnetic field 
is set at the level crossing and the microwave frequency is set at the center of the 
/5-1/2 resonance we have [30]: 

(da \ ^ 

Two quantities are useful in assessing the sensitivity of this type of APV experiment. 
The ffist is the fraction of the rate that is due to the PNC interference term, which we 
will call the asymmetry A: 

A = A^^Cro (18) 
72p Atoba. 

-1.1 X 10"^ , , 

tan 

where we used C\r, = —0.459, and b^/ex = tan0. The asymmetry can be set 
experimentally by choosing the angle 0. Here we just assume the value tan0 = 0.1 
which gives an asymmetry of —1.1 x 10~®. 

A more important quantity for assessing the viability of H/D APV experiments 
is the counting time required under ideal conditions, i.e. assuming one can achieve the 
shot-noise limit, and that the transition rate is well above background. The signal count 
5* is proportional to {Ap^f^Apc + c.c.) times the factor rjJT, where rj is the detection 
efficiency for the /5_i/2 state, J is the metastable beam intensity per hyperfine component 
(atoms/s), and T is the total counting time. The noise count is proportional to the 
square root of ?7JT|Apcp. Using Eqs. [131 [161 and [T7I we find: 

— = 4 V^CiD J TYT^ ivJT) I . 20 
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Note that Eqs. [2D] does not depend on and, in the approximation we use here, the 
counting time is independent of the asymmetry chosen for the experiment. In fact, as 
pointed out by Hinds[5l], if the phase between Apc and ApNc is chosen to maximize 
the interference term, the signal-to-noise is independent of Apc- 

Equation [20] has a maximum with respect 'ja_f_i/2 given by the condition: [54j : 

7a+,/, = 2.5/r, (21) 

which, for a given r, can be satisfied by adjusting the microwave power. Since the 
required power is readily obtained for the experiments considered here, we use Eq. [21] 
to eliminate la+^/i from Eq. [20] with the result: 

(J^^ =1.8\V^C^d\{vJTt/^2p)'^'. (22) 

The definitions of the parameters in this equation are summarized in Table [3j Eq. [22] 
indicates that, apart from maximizing the beam intensity and detection efficiency, 
making r as large as possible is the key requirement for maximizing S/N. This means an 
experiment based on a slow beam and a long interaction region. A hidden assumption 
here is that the linewidth of the a+1/2 — > P-1/2 transition is of order the width of the 
a+1/2 state (Eq.[T7|). In particular, the magnetic field must be constant and homogeneous 
at the level: 

M< -'-^-'^ ^l^. (23) 

B t^a+i/2/3_i/2 ^ 

Eq. [22] also applies to measurements performed at the other level crossings listed in 
Table [21 providing the combinations of coupling constants listed there are substituted 
in place of Cir,. We can use the optimized signal-to-noise ratio to estimate the 
minimum limits that can be placed on the electron-nucleon, weak-neutral-current 
coupling constants from experiments in H and D. To do this, we assume that the counting 
time T should not be more than about 5 days; with the understanding that, in reality, 
considerably more counting time would be required because of backgrounds and the need 
to study systematic errors. Then, taking the overall detection efficiency to be 77 = 0.5, 
we plot the percent errors in the four constants as a function of the metastable beam 
intensity per hyperfine component J, and the time of passage of an atom through the 
interaction region r. We also give the magnetic field homogeneity required. The results 
for the deuterium couplings are given in Fig. [HI and those for hydrogen are given in Fig. [7] 
For example. Fig. [6] indicates that to measure Cid to ~ 0.3% with an interaction time 
T ~ 100/is, requires a beam intensity of at least J ^ 5 x lO^^s"^ and AB/B ~ 10~^. 
For longer interaction times r, the requirement on the beam intensity can be reduced 
but this also leads to more stringent specifications on the B-field homogeneity. The 
same experimental conditions would provide a measurement of C2D to better than 20% 
and (based on Fig. [7]) a ~ 4% measurement of Cip. "Interesting" experiments can be 
done with smaller beam fiuxes as can be seen from the lines in Figs. [6] and [7] indicating 
the model-independent limits on the four coupling constants as well as the "dot-dashed" 
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Figure 6. Deuterium APV sensitivity. Solid lines indicate projected minimum 
experimental uncertainties in Cid (in percent of the SM value) for a generic APV 
experiment in D based on Eq. [22l These are plotted as a function of the time an atom 
spends in the interaction region r, and the beam intensity per hyperfine component 
J. The values assumed for the other parameters in Eq. [221 are summarized in Table [H 
The dashed lines are the projected uncertainties for a C2D measurement at the /3e 
crossing. The dot-dash line represents the uncertainty needed in a measurement of 
CiD to determine sin^ 6^ with a precision equal to the current APV limit of 0.9%. The 
lines marked (7mi(C2_D) and crmi(C'i£)) represent the uncertainties in these parameters 
based on a model-independent-fit to all neutral current data|T|. The axis on the right 
gives the magnetic field homogeneity AB/B based on Eq. [23l 




Figure 7. Hydrogen APV sensitivity. Projected minimum experimental uncertainties 
(as a percentage of the SM values) in the hydrogen weak interaction coupling constants 
based on Eq. [22l and plotted as a function of r and J. (See caption to Fig. [6l) The 
solid lines refer to a measurement of Cip and the dashed lines refer to measurement 
of C2p at the /3e crossing. The dot-dashed line corresponds to the uncertainty in 
Cip required to match the current APV limit on sin^ 9w The lines marked <Tini(Cip) 
and (Tmi(C'2p) correspond to the uncertainties in these parameters determined by the 
model-independent-fit discussed in Ref. [T]. 
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lines which show the experimental uncertainty required to match the current APV limit 
on determination of sin^ 9^. 

6. Systematic Effects 

The 2^Si/2 and 2^Pi/2 levels in hydrogen can be mixed by unwanted electric fields and 
this can mimic parity mixing. The closeness of the levels enhances the PNC effects but 
it also increases the sensitivity of the experiments to systematic effects due to unwanted 
electric fields. In this section, we will consider two of the most troublesome effects 
expected in the experiment discussed in Section [51 

A number of essential elements of the treatment of systematic errors will not 
be covered here, however. These include the techniques for accurately aligning the 
apparatus, and reversing the vectors in the interaction region including: the polarization 
0, the microwave propagation k, and the magnetic field B. We will assume that 
the apparatus can be aligned to the levels achieved in the former experiments e.g. 



However, in [Appendix C| a scheme for continuously monitoring the microwave 
polarization is discussed, that could make this reversal exact to the level of statistical 
accuracy. Another crucial issue concerns effects caused by variations in the intensity or 
direction of the fields as the atoms move through the apparatus |18]. These effects will 
not be discussed here. They are addressed in detail by e.g. Levi and Williams [H]. 

To understand the problem of unwanted electric fields, we introduce electric fields 
i?", Ey and E^ into our amplitudes. We also allow a misalignment of the microwave 
field which gives components and b^ in addition to the desired fields. The magnetic 
field defines the z axis and we are free to take the microwave propagation vector k to 
lie in the y — z plane. With these definitions, after adding the most important electric- 
field- induced terms, A becomes |55] 

2hA = -(b^ + ib%)/io cose - i(e^ + ie'^)dwA(^_,^^f_^^, 

+ {E: + \El)eA^p.y,e^y, - E^e. + ie^)dl^p_,^^f_,^,. (24) 

Here, dw = VCiDd2. Under a polarization reversal — » —(f), the microwave components 
ez, e^, and b^^ change sign while e^, b™, and b^ do not. 

The first two terms in Eq. [2l] (with = = 0) are the amplitudes Apc and Ap^c 
discussed in the previous section and the last two terms are the electric field induced 
amplitudes which are of concern in this section. The first point is that interferences 
among the electric field induced terms {ElEl interferences) are relatively small since 
there are no applied electric fields. Also, these terms are readily eliminated because they 
are even in the reversal of k, whereas the ElMl PNC pseudoscalar term is odd under 
this reversal. Of the i^lMl interferences, those that change sign under a reversal of 
the polarization are of most concern. The two most significant effects are: (1) a term 
caused by a transverse electric field which we will call the motional field effect, and, (2) 
a term caused by an electric field along B which we will call the stray E^ systematic 
effect. In the following we estimate the size of these effects. Based on Eq. [2T]and Eq. [171 
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and taking r = 300/iS and tan 0=0.1, we find: 

^ 0.16V /cm, e, 0.016V /cm, ^ 5 x lO'^G. (25) 

6.1. Motional electric field systematic effect 

The third term in Eq. [2l] arises from mixing of /5_i/2 with e+1/2 by a transverse electric 
field and, unlike the PNC term (2"'^ term in Eq. [2^ it is not resonant at the f3f crossing. 
The largest term is proportional to the combination of fields, b™i?^. This involves 
an electric field transverse to the magnetic field. The most significant field of this kind 
arises from the v X B "motional" electric field. 

It is of interest to form the ratio of the interference term arising from the 
motional electric field to that of the pseudoscalar interference term. This is: 

j^mot ^ 8 ^ ^Q5 (-g"';^)sinx ^ ^26) 

where x is the microwave misalignment angle defined by b™ = b^; sin x, and {E"^"^) is 
the motional field (in Volts/cm) averaged over the beam velocity distribution. If the 
transverse velocity distribution were cylindrically symmetric and the axis of symmetry 
exactly coincident with a homogeneous magnetic field, the motional field systematic 
would vanish. For simplicity, we represent the imperfections of the real apparatus by a 
misalignment between the beam velocity and the magnetic field. Based on Fehrenbach's 
experiment [44] we take this angle to be 0.1 mrad. For a slow (77K) deuterium beam, 
and a magnetic field of ~ 1200 G, we find {E^°^) ~ 10~^^/cm. Assuming the microwave 
misalignment angle sinx = 10~^, we have: 

^mot ~ Ml ~ 0.02 (27) 

The motional field effect can be distinguished from the pseudoscalar term because it is 
sensitive to slight changes in the direction of the B-field and it is non-resonant at the 
crossing. 

For a new deuterium APV experiment, it would be desirable to improve the 
symmetry and alignment of the beam over what has been achieved in the past in order to 
suppress this effect by another order of magnitude so that it is smaller than 0.3% of the 
pseudoscalar term. But, regardless of the level of suppression obtained, this effect can be 
accurately measured and a correction can be applied to the final result of the experiment. 
To do this, the experiment must incorporate a continuous monitoring of the motional 
electric fields while the APV data are being taken. The general method for doing this 
is well known, as it is an essential technique in APV experiments. The unwanted fields 
are monitored by using the atoms themselves [371 EG]. Iii the present case, the motional 
field could be monitored by, e.g. periodically moving the magnetic field off the center 
of the level crossing (retuning the microwave frequency to an adjacent resonator mode). 
Away from the crossing, ApNc is suppressed, so the signal is dominated by the motional 
electric field effect. The limit on the motional field effect is determined by the statistical 
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accuracy of the on-line monitoring, which by the arguments of Section \5\, again reduces 
to the need for an intense, slow metastable beam. 

The largest systematic error in the measurement of Fehrenbach was caused by a 
motional electric field and it is natural to ask why this is less problematic here. The 
first point is that the velocity of the slow D beam is a factor of ~ 300 smaller than 
his fast beam so the motional field is smaller by the same factor. Also, his experiment 
was performed at the crossing of the e and f3 levels which are mixed by the (transverse) 
motional field, so for an experiment done at the (3f crossing there is an additional 
relative suppression of approximately 10. The net result is a suppression of a factor of 
approximately 3000 for the motional field effects in the experiment discussed in Section |5l 
relative to those in Fehrenbach's experiment. 



6.2. Stray systematic effect 

The last term in Eq. [21] arises from mixing between the (3-i/2 and f-1/2 states by an 
unwanted "stray" electric field parallel to the magnetic field. Such a field could arise 
from surfaces in the interaction region, charged particles in the beam, etc. We call the 
most significant term, proportional to the fields h^E'^, the stray E^ systematic effect. 
The microwave field combination is the same as that for the pseudoscalar interference, 
but this term has a dispersive shape as a function of magnetic field. It goes to zero near 
the center of the P-i/2f-i/2 crossing, and has a maximum at i^/3_i/2/_i/2 ~ 72p/2. The 
exact magnetic field where the stray Ez term vanishes depends on couplings with more 
distant levels and it can be calculated to high precision. Due to field inhomogeneities, 
frequency shifts, the finite linewidth of the transition, etc., one cannot completely zero 
the stray E^ effect. Rather, we introduce a parameter 5ab which is the degree of 
suppression of this term from its maximum value obtained by setting the magnetic field 
to the proper value. There is an additional suppression of this term since it is odd under 
a reversal of B, while the pseudoscalar is even in this reversal. We assign another factor 
to account for this suppression. 
The ratio R^'' of the stray E^ interference term to the pseudoscalar is found using 
the same microwave field strengths as previously: 

CiD 

[El) is the average electric field along z seen by the atoms, expressed in V/cm. 

The field {Ef) can be minimized by careful design of the interaction region, but 
here we just assume the value achieved by Fehrenbach which was {E^) ~ ImV/cm. We 
take the suppression factors to be = ^ab = 10~^, based on his results. Using these 
values we find: 

1 n-3 

R^^ ^ —— ^ -2 X 10"^ (29) 

where the coupling constant was taken from Table [2J This estimate is an order of 
magnitude smaller than the motional field effect, and is smaller than our goal of 0.3% 
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uncertainty in a measurement of the pseudoscalar term. The stray systematic effect 
is readily measured while the APV data are being recorded by setting the magnetic 
field off the crossing and looking for a term which is polarization dependent, odd in 
the magnetic field, and switches sign when the field is tuned to the other side of the 
crossing. 



7. Summary 

The major argument for consideration of parity experiments in H and D, that there 
is negligible atomic physics uncertainty in the interpretation of the results, remains 
compelling. Parity violation experiments in deuterium could play a role in looking 
for new physics by providing a low energy measurement of sin^^vK to confirm the SM 
prediction for the "running of sin^6'iy"- A 0.3% measurement of Cm would improve the 
limit on sin^^iy provided by APV experiments by a factor of three. 

Hydrogen and deuterium APV experiments could also contribute to determination 
of the electron-nucleon weak neutral current coupling constants listed in Table [H Only 
CiD is accurately constrained by current data. The uncertainty in Cip is about the size 
of the SM prediction, and C2p and C2D are completely undetermined. The values of the 
individual constants could be determined by measuring APV at several level crossings. 
An updated list of the SM predictions for the parity violating matrix elements at level 
crossings in H and D are given in Table El 

In Section [5] we analyzed a generic deuterium APV experiment, and found that a 
0.3% measurement of Cid would be possible if a slow metastable beam with an intensity 
of ~ 5x10^^ deuterium atoms/s per hyperfine level were available. In Appendix B we 
discuss some ideas for obtaining such a beam. 

In Section [6], we discussed two systematic effects and found that, although they 
provide a challenge for a deuterium APV experiment, they are not insurmountable 
obstacles. After designing the experiment to minimize the systematic effects due to 
unwanted electric fields as much as possible, the elimination of these effects is done by 
monitoring them during the APV data taking, using the atoms themselves to measure 
the fields. Thus, the limitation on the experimental sensitivity due these systematic 
effects is dependent on counting statistics, and, in this sense, also depends on the 
intensity of the metastable beam. 
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Appendix A. Metastable quenching from beam-beam collisions 

Metastable quenching from beam-beam collisions places a practical limit on the density 
of metastable D(2s) atoms in a beam suitable for a deuterium APV experiment. Forrey 
et al. [57j calculated the cross sections for ionization, excitation-transfer, and elastic 
scattering in collisions between metastable hydrogen 2s atoms at thermal energies. The 
most important metastable destruction mechanism was found to be double excitation 
transfer to H(2p) : 

H{2s) + H{2s) ^ H{2p) + H{2p) (A.l) 

for which they found a cross section a'^^ = 9 x lO^^^cm^ at = 4.1 meV, varying as 
E^^^"^ at higher energies. 

If we require that the mean free path A^/p of the D(2s) atoms be a meter, this puts 
an upper limit on the metastable density: 

P2S < ~ (A.2) 



mfp 



Then, for a beam with an average velocity of 10^ cm/s, the D(2s) flux density must be 
< 10^^ cm~^s~^. We can take this limit to apply to an individual hyperflne component 
since one can quench out all unwanted components immediately after formation of the 
metastable beam using a deuterium spin fllter[58]. For a metastable beam of SxlO^'^s"^ 
per hyperflne component, the beam diameter has to be > 2.5 cm. 

The metastable beam can also be quenched from beam-beam collisions with D(ls) 
and D2. The cross section[59] is about 2 x lO^^^cm^, so for a mean free path of 1 m 
at 10^ cm/s, the flux density must be < 5 x lO^^cm^^s^^. For a 2.5 cm diameter beam 
this corresponds to 2 x 10^'' D(ls)/s. 

Appendix B. Optical production of metastable hydrogen and deuterium 
beams 

One possibility for obtaining a slow metastable D(2s) beam with suflicient intensity 
for a deuterium APV experiment is to use optical pumping starting with an intense 
dissociated deuterium beam. One such beam was produced by Harvey [60]. In his 
scheme, a beam of hydrogen atoms in the ground state was pumped to 3p using a 
Lyman-/? discharge lamp. The 3p level decays to the 2s state about 12% of the time. 
This work produced a beam with a flux of only 2 x 10^ metastable atoms/s because 
of the low intensity of Lyman-/5 radiation from the lamp. The development of intense 
beams of UV radiation using FEL technology might provide a considerable improvement 
in this respect. 

The photon absorption cross section [61] for the Is 3p assuming the radiation 
does not resolve the 3p flne structure is given by 

3A^73pis , , 

= -7^ (B.l) 

27r 7iine 
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where A is the transition wavelength, 'jspu is the ground state decay rate for the 3p level 
and 7iine is the transition linewidth. We find: 

(^is3p ~ 7— --cm^ B.2 

We are interested in the amount of UV power at 102.6 nm needed to produce a 
metastable D(2s) beam (in a single hyperfine component) of 5 xlO^^ D(2s) s~^ to serve 
as a benchmark for FEL development. Assuming a ground state beam of 2 x lO^^s"^ 
with a diameter of 2.5 cm and a UV beam with a 3 cm diameter and accounting for 
the 12% branching ratio for 3p 2s, we find the power per bandwidth required is 40 
mW/GHz. 

If intense Lyman-a radiation is available, metastable hydrogen can be obtained by a 
three step process; two photon pumping Is ^ 2p ^ 4d followed by decay to 2s. Lyman- 
a radiation drives the transition Is —>■ 2p which is followed by a 2p — > 4d transition 
saturated by a visible laser. Atoms in 4d have a 3% branch to the 2s state |62J. This 
scheme is similar to that used by Young, et al. for production of metastable Kr atoms 
[63t [M] . Although the branching ratio for the 4d state is not as favorable as for the 3d 
state, the oscillator strength for the transition Is — > 2p {fab = 0.416) is five times larger 
than for the Is 3p transition. Also, an FEL capable of producing several watts of 
Lyman-a radiation is under development at Jefferson Lab in Virginia [65j. 

Another possibility for production of a metastable H beam is based on the two- 
photon transition Is — > 2s with 243 nm radiation [661 EH EH [69] . This is a weak process 
since it does not involve an intermediate resonance state, and requires high intensity 
radiation. The needed intensity is normally obtained by focussing the light to a small 
spot, but this is incompatible with the limitation on the density of the metastable 
hydrogen atoms discussed in Appendix A which requires a photon spot size of at least 
5 cm^. For this reason, two-photon excitation does not appear to be a viable option 
for production of the required metastable beam. This conclusion is also supported by 
existing experimental data [70l [71] . 



Appendix C. Normalization Transition 

In the generic experiment discussed in Section [5], the psuedoscalar interference term is 
measured by picking out the part of the a+1/2 — > (^-1/2 transition rate that changes sign 
under the polarization reversal (f) —>■ —(f). Ideally, the microwave field that drives the 
PC transition would have the same strength in each polarization state, so the PC part of 
the transition would average to zero over the course of the experiment. In practice, it is 
not possible to achieve this to sufficient accuracy, and several independent reversals are 
used[41j. However, it is possible to monitor the field strength in each polarization state 
during accumulation of APV data, and make this reversal exact within the statistical 
uncertainty of the monitor. One way to do this is by periodically switching on an 
auxiliary static electric field E'^ perpendicular to B. This adds a term to Eq. [13] of the 
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form: 

2hA^ = rf^A,_,/,e,,/,(i?f ) + ii?De. (C.l) 
The field direction can be chosen so that the new term is nearly orthogonal to the PNC 
amplitude. A^r is proportional to the microwave electric field component which is 
directly proportional to the microwave magnetic field component b^^ that needs to be 
measured. (See Fig. 1). The field strength E^^^ is chosen so that A^v is much larger 
than A (Eq. [T3ll so that, when E^'^^ is on, Ajy dominates the a+1/2 P-1/2 transition 
rate, and the contribution of the PNC interference term is sufficiently small that it does 
not contribute to the measurement of the field strength. 
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